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We describe the process of parametric amplification in a directional coupler of quadratically 
nonlinear and lossy waveguides, which belong to a class of optical systems with spatial parity-time 
(PT) symmetry in the linear regime. We identify a distinct spectral parity-time anti-symmetry 
associated with optical parametric interactions, and show that pump-controlled symmetry breaking 
can facilitate spectrally selective mode amplification in analogy with PT lasers. We also establish 
a connection between breaking of spectral and spatial mode symmetries, revealing the potential to 
implement unconventional regimes of spatial light switching through ultrafast control of PT breaking 
by pump pulses. 
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Light propagation in waveguiding structures with spa¬ 
tially distributed sections of loss and gain can be anal¬ 
ogous to quantum wavepacket dynamics governed by 
a parity-time (PT) symmetric Hamiltonian [T]. Be¬ 
low a certain gain/loss level, such systems support PT- 
symmetric optical modes, which then exhibit the same 
average loss or gain laia. However when gain or loss 
is increased, the PT-symmetry of modes breaks, and a 
mode with the strongest gain (or smallest loss) domi¬ 
nates, as demonstrated experimentally dill- The phase 
transition associated with such PT-symmetry breaking 
opens new possibilities for light manipulation, such as 
PT-symmetric lasers [3 E]. Such lasers can achieve 
single-mode operation, where small difference in medium 
gain leads to a dramatic difference in mode amplification 
below and above the PT breaking threshold. 

Parametric amplifiers are commonly used as an inte¬ 
gral part of optical setups enabling flexible wavelength 
conversion and tunable signal gain, extending the range 
of lasers where gain media are limited to particular wave¬ 
lengths. Wave amplification is efficiently realized in the 
regime of difference-frequency generation in media with 
quadratic optical nonlinearity [7]. Importantly, the am¬ 
plification rate is determined by the pump, enabling ul¬ 
trafast all-optical tunability. In this work, we reveal the 
potential of PT-symmetric systems for optical parametric 
amplification, and identify a new regime of spectral PT 
anti-symmetry in such devices. Such devices can, on one 
hand, realize ultrafast spatial signal switching through 
pump-controlled breaking of PT symmetry, and on the 
other hand enable spectrally-selective mode amplification 
in analogy with PT lasers. 

We consider a directional coupler composed of two 
waveguides in quadratically nonlinear medium, where 
modes exhibit different loss in each waveguide. It can be 
realized experimentally based on LiNbOs couplers where 
parametric gain was demonstrated [8], as well as other 
platforms with nonlinearity. The loss can be intro¬ 
duced, for example, by depositing a thin layer of metal on 
top of th waveguide [1] . An illustration of such structure 


with loss in one waveguide is presented in Fig. In the 
linear regime, at low light intensities, such coupler real¬ 
izes PT-symmetric optical system [1]. However at higher 
intensities the effect of quadratic nonlinear interactions 
becomes important. It was predicted that in the regime 
of second-harmonic generation, when the signal and idler 
waves have identical spectra at half of the pump fre¬ 
quency, parametric interactions in PT couplers can sup¬ 
port a rich variety of nonlinear modes [9]. Furthermore, 
the formation of quadratic solitons in spatially extended 
PT-symmetric structures was analyzed in detail mm- 
However, the fundamentally important regime of para¬ 
metric amplification in PT systems remained unexplored. 

We analyze the process of optical parametric amplifi¬ 
cation based on nonlinear mixing between strong pump, 
signal and idler waves, as illustrated in Fig.[2 We model 
the wave propagation using coupled-mode equations [7] 
in the regime of narrowband and undepleted pump, 
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Here the subscripts stand for signal (‘s’), idler (‘i’), and 
pump (‘p’) waves in two waveguides (‘1’ and ‘2’), E are 
the mode amplitudes, ^ is the propagation distance along 
the waveguides, (3 are the propagation constants (for the 
pump mode, the propagation constant is adjusted to ac¬ 
count for quasi-phase-matching through periodic poling 
of the ferroelectric domains 0 ), 7 are the linear loss co¬ 
efficients, C are the mode coupling coefficients between 
the waveguides, and x are effective quadratic nonlinear 
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FIG. 1 . Scheme of PT-symmetric nonlinear coupler with 
linear absorption in one waveguide. 


coefficients. 

We find that PT-symmetry in the regime of paramet¬ 
ric amplification can be achieved in the near-degenerate 
case when the waveguide parameters are practically the 
same at signal and idler frequencies, i.e. 7^1 = 7^1 = 71 , 
7 i 2 = ls 2 = 72 , and Cg = Ci = C. We assume that 
the waveguides are engineered such that ^si = l3s2 = 
and Pii = Pi 2 = A, which is the regime required for 
linear PT-symmetry mill. Under usual experimental 
conditions the mode at higher pump frequency is much 
stronger localized compared to signal an idler m, lead¬ 
ing to suppressed coupling between the waveguides and 
also very small sensitivity to metal deposited on top of 
waveguides, meaning that Cp = 0 and 7^1 = 7^2 = 0 . 
We further consider a case of equal pump propagation 
constants in two waveguides, f3pi = Pp 2 = Pp. The latter 
condition does not need to be satisfied if pump is cou¬ 
pled to one waveguide, which already enables full range of 
mode switching and amplification control as we demon¬ 
strate below. 

To reveal the PT-symmetry properties of Eqs. Q, we 
represent the equations for signal and idler waves in the 
Hamiltonian form, 

i- = Ha, ( 2 ) 
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FIG. 2. (a,b) The largest mode gain (white line marks 

a zero gain level) and (c,d) number of PT-symmetric mode 
pairs vs. the input pump in the first waveguide and the 
phase-mismatch. The linear losses are (a,c) 72s = 72z = 1 , 
(b,d) 72s = 72 z = 3 . For all the plots 71s = 71^ = 0 , C = 1 
and A2 = 0. 


Here T is a time-reversal operator which changes z ^ — z 
and performs a complex conjugation. The parity opera¬ 
tors operate in speetral domain^ interchanging the signal 
and idler waves. 


Vi^± — {Ogi , 7 ^ 2 ,± — {cis2 ^ • (^) 

We define the parity operators with both symmetric (‘+’) 
and antisymmetric (‘-’) transformations, since the latter 
will be useful in the following analysis. We note that 
such unusual symmetry is fundamentally different from 
the previously studied antisymmetric PT-metamaterials 
with modulated dielectric and magnetic properties m- 
Since the Hamiltonian is linear in the undepleted pump 
regime, the dynamics of signal and idler waves is defined 
by the eigenmode solutions, 

a(2;) = a(cr) exp(icr2;), ( 7 ) 
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where P = {Pp — Ps — Pi) j2 defines the phase mismatch 
of parametric wave mixing, and Hi ^2 are the normalized 
input pump amplitudes. Hi = — 0 ) and H 2 = 

X2Ep2{z = 0 ). 

A key result of our analysis is that the Hamiltonian 
possesses a speetral anti-PT symmetry^ corresponding to 
a negative sign on the right-hand side of the following 
equality, 

(5) 


where a is an eigenvalue. The real part. Re (a), defines 
the phase velocity, whereas the imaginary part deter¬ 
mines the modal gain, F = —Im(cr). 

We determine the effect of spectral anti-PT symme¬ 
try on the eigenmode properties by substituting Eq. 0 
into Eq. 0 and applying PT operator. We obtain that 
if a(cr) is an eigenmode, then a(—a*) = 7^i^+7^2,+Ta(cr) 
is also an eigenmode. There are two possibilities how 
these relations can be satisfied. Eirst, the mode can be 
PT-symmetric, when PT transformation maps the mode 
profile to itself (up to an overall phase), which happens 
if —cr* = cr, and accordingly Re (a) = 0. Such spee- 
trally PT-symmetrie modes would generally experienee 
gain/loss different from other modes^ since there are no 


Vi,+V2,+Tn = -nvi,+V2,+T. 
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FIG. 3 . Mode eigenvalues vs. the pump amplitude in the 
first waveguide (A2 = 0): (a) negative imaginary part defining 
mode amplification, F = —Im(cr), and (b) real part defining 
propagation constant, Re(cr). Parameters are 71s = = 0 , 

72s = 72z = 1, C' = 1, and /3 = 0. 


specific relations for the value of Im(cr). Second, if the 
mode profile has broken PT-symmetry, the PT transfor¬ 
mation relates two different modes with eigenvalues a2 = 
—( 7 \. It follows that a pair of spectrally PT-broken modes 
experience the same loss or gain, Im(cri) = Im((72), but 
they have different phase velocities, Re(cri) = —Re(cr2). 
Remarkably, the established relations of mode symme¬ 
try and gain/loss are reversed in comparison to previ¬ 
ously studied spatial PT-symmetry in directional cou¬ 
plers mm, due to the spectral anti-PT symmetry of 
parametric wave mixing. 

We now demonstrate that the modal PT-breaking can 
be controlled by the pump beam. Due to the elec¬ 
tronic nature of quadratic nonlinearity, such tuning can 
be ultrafast, directly following the pump profile in real 
time. Overall, there are four eigenmodes of the Hamilto¬ 
nian. Accordingly, there can be three possible symmetry 
regimes: (i) there are two mode pairs with broken PT- 
symmetry, (ii) one pair of PT-broken modes and a pair of 
PT-symmetric modes, or (iii) two pairs of PT-symmetric 
modes. As an example, we present numerical analysis 
of the mode properties in a coupler with loss in the sec¬ 
ond waveguide, and pump coupled to the first waveg¬ 
uide. We show the largest mode gain F = Max[—Im(cr)] 
in Figs. [^a,b) and the number of PT-symmetric mode 
pairs in Figs. |^c,d), depending on the pump amplitude 
Ai and the phase-mismatch, for different values of lin¬ 
ear loss. We observe that in the regime when all modes 
have broken spectral PT symmetry [blue shaded regions 
in Figs. [ 2 ];c ,d)], the modes experience negative gain [c.f. 
Figs. [^a,b)]. This happens because pairs of eigenmodes 
exhibit the same amount of gain/loss, and effectively the 
amounts of gain and loss are averaged out between the 
eigenmodes. However upon transition to the region with 
spectrally PT-symmetric modes [green and red shaded 
regions in Figs. [^c,d)], there appears an unequal redis¬ 
tribution of gain and loss between the modes. One PT- 
symmetric eigenmode exhibits gain much larger then all 


other modes, while the latter experience stronger loss. 
Such sensitivity of amplification to PT-breaking thresh¬ 
old could be used to discriminate between multiple spec¬ 
tral modes, analogous to the concept of PT-lasers 016 ]. 

Next, we establish a connection between the spectral 
PT-symmetry identified above, and the spatial wave dy¬ 
namics due to waveguide coupling. We first consider 
a special case which can be treated analytically, corre¬ 
sponding to perfect phase-matching, /3 = 0 , and addi¬ 
tionally Im(A^A2) = 0 . The latter condition can be 
transformed to Im(Ai) = Im(A2) = 0 under a substi¬ 
tution a ^ aexp(i(/:?) with appropriately chosen constant 
phase ip. Then, we find that 

Vi,+V2,-n = nvi,^V2,-, (s) 

Vi^-V 2 ,+n = nvi^-V 2 ,+. ( 9 ) 

Accordingly, there appear two pairs of eigenmodes with 
= rjOsi and a*2 = where one pair with 

?7 = +1 has a profile symmetric with respect to Pi^+p2,- 
and second pair with 77 = —1 conforms to the symme¬ 
try Vi-V2,+> Remarkably, the signal dynamics of such 
modes is governed by equations resembling those for a 
linear PT-symmetric coupler mm, 
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( 10 ) 

where a^ = {asi;as2)- We see that the linear loss is mod¬ 
ified due to the parametric gain determined by the pump 
amplitudes. The coupler Eqs. ( 10 ) are symmetric with 
respect to spatial PT-symmetry, up to a gauge transfor¬ 
mation [1] expressed through the identity matrix (X) in 
the following relation. 


Vl2r{nr -pT) = {Hr -pT)Vl2T, (11) 


where the spatial parity operator 1^12 swaps the mode am¬ 
plitudes between the waveguides, a si ^ as2 and an ^ 
ai2. The coefficient p = (77A1—71)/ 2 + (—77^2—72 )/2 de¬ 
fines the average gain or loss between the two waveguides, 
which depends on the pump amplitudes and signal/idler 
mode symmetries corresponding to the different signs of 
‘77’. The mode eigenvalues are 

a =-ip ± (1/2) 7^2 - [71 - 72 - r/(Ai + ^2)]"- (12) 

We present characteristic dependencies of the eigenval¬ 
ues on the pump amplitude in the first waveguide in 
Fig. i We find that both spatial and spectral PT- 
symmetry breaking occurs simultaneously at the threshold 
I71 — I2 — v{Ai + ^2)1 = C. However, the spatial and 
spectral symmetries are opposite: a mode pair is spatially 
PT-symmetric and has spectrally broken symmetry be¬ 
low threshold, whereas the situation is reversed above 
the threshold. We show the evolution of the signal inten¬ 
sity along the waveguides is shown in Fig. At lower 
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FIG. 4. Evolution of the signal mode intensity along the 
waveguides: solid (green) lines — the first waveguide with¬ 
out loss, dashed (blue) line — the second waveguide with 
loss. The pump amplitude is (a) below PT-symmetry break¬ 
ing threshold {Ai = 0.5) and (b) above threshold {Ai = 1.5). 
Parameters are 71 s = ju = 0, 72 s = 72 z = 1, C = 1, and 
13 = 0 . 


pump powers [Fig. [^a)] signal periodically switched be¬ 
tween the waveguides due to the beating between two 
modes exhibiting the same (negative) gain. However for 
stronger pump above the PT threshold [Figs.|^b)], only 
one supermode with the strongest gain dominates and 
accordingly oscillations are suppressed. 

Finally, we investigate numerically a connection be¬ 
tween spectral PT-breaking and spatial mode dynam¬ 
ics in a general case of non-zero phase mismatch. We 
present in Fig.j^the fraction of signal intensity in the first 
waveguide depending on the phase-mismatch and the 
pump amplitude for different propagating distances. We 
find that spatial dynamics strongly depends on the spa¬ 
tial PT-symmetry of modes in the linear regime (at low 
pump amplitudes Aj ^ 0 ). If the linear modes are PT- 
symmetric, I71 — 72I < C, then increase of pump ampli¬ 
tude can control the period of mode coupling between the 
waveguides, while the oscillations get suppressed close to 
the spectral PT threshold, see Figs.|^a,c) and Fig.j^c). 
However if the linear modes are PT-broken, I71 — 72 1 > C, 
then mode beating between the waveguides still occurs 
below the spectral PT threshold, but with very small 
modulation amplitude, see Figs. [^b,d) and Fig. |^d). In 
all cases, the mode beating changes faster with pump 
amplitude at longer distances. 

In conclusion, we have identified an anti-PT spectral 
symmetry of a parametric amplifier based on a quadrati- 
cally nonlinear coupler with different losses in two waveg¬ 
uides. For pump powers below the threshold, the modes 
form pairs with broken PT symmetry and same gain/loss, 
whereas above the threshold one PT-symmetric mode ex¬ 
periences the largest gain. This can facilitate spectrally- 
selective mode amplification, and we expect that this 
effect will become even more pronounced in resonator 
structures in analogy with single-mode operation of PT 
microring lasers. We have further established an under- 



FIG. 5. Fraction of signal intensity in the first waveg¬ 
uide, |asi|^/(|asi|^ + |as 2 |^), vs. the input pump in the first 
waveguide and the phase-mismatch. Propagation distances 
are (a,b) L = 5C and (c,d) L = IOC. Parameters correspond 
to (a,c) Figs.|^a,c) and (b,d) Figs. | 2 i:b,d). 


lying connection between spectral anti-PT and conven¬ 
tional spatial PT mode symmetries, which reveals the 
possibility to control spatial light switching and ampli¬ 
fication through parametric gain. Accordingly, the sug¬ 
gested platform can implement various unconventional 
regimes of light control previously suggested for PT- 
symmetric structures, such as unidirectional m and 
nonreciprocal m operation, but with the advantage of 
ultrafast all-optical control of PT-symmetry regimes by 
pump pulses. We anticipate that, due to the univer¬ 
sality of parametric amplification processes, these con¬ 
cepts will be extended to different physical mechanisms 
including Kerr-type optical nonlinearity, as well as a 
broad range of other systems including cold atomic Bose- 
Einstien condensates m in engineered potentials [num 
and exciton-polariton condensates [ 19 ]. 
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